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Abstract 

The behaviour of a particle with a spin 1/2 and a dipole magnetic 
moment in a time- varying magnetic field in the form (hocn{ujt, k), 
hQsn{u;t, k) , HQdn{ujt, k)) , where lo is the driving field frequency, t is 
the time, /iq and Hq are the field amplitudes, cn, sn, dn are Jacobi 
elliptic functions, k is the modulus of the elliptic functions has been 
considered. The variation parameter k from zero to 1 gives rise to a 
wide set of functions from trigonometric shapes to exponential pulse 
shapes modulating the field. The problem was reduced to the solution 
of general Heun' equation. The exact solution of the wave function 
was found at resonance for any k. It has been shown that the transi- 
tion probability in this case does not depend on k. The present study 
may be useful for analysis interference experiments, improving mag- 
netic spectrometers and the field of quantum computing. 
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1 Introduction 



Rabi llj studied the temporal dynamics of a particle featuring a dipole 
magnetic moment | in a constant magnetic field Hq, directed along the z- 
axis, and another varying magnetic field = hocosut, Hy = hosinut ro- 
tating with a frequency u perpendicular to Hq { Hq, ho are the amplitudes). 
There are several methods of modulating magnetic fields while studying the 
phenomenon of magnetic resonance [2 . This work focuses on the temporal 
evolution of a particle with a dipole magnetic moment in a distorted mag- 
netic field described by H{t) = {hocn{ujt, k), hosn^ut, k), Hodn{ut, k)). Such 
field modulation under a changing modulus k of the elliptic functions from 
zero to unity describes an entire class of field shapes from trigonometric P] 
to pulsed exponential OEj. 

The Schrodinger equation of a wave function that describes the dy- 
namics of a particle featuring a spin of | and a magnetic moment in a time- 
varying magnetic field H{t) is given 

zhd,^{t) = ^aHm{t), (1) 

where g is the Lande factor, fiQ is the Bohr magneton, and the Pauli matrices 
are 

a = {o-xiCTyjCTz) = , ^ , ^ ° . The solution to the 

Schrodinger equation could be found by expanding in eigenf unctions of ma- 
trix az 

^ (t) = Q (t) + (^^ (t) . (2) 

Functions \E'i {t) and ^2 if) describe the probability amplitudes without a 
spin flip and with one, respectively, and by definition obey the normalization 
restriction 

i^i(t)|2+i^2(t)r = i (3) 

and are assumed to be defined at the initial time t. 
Considering the condition (2), Eq. (1) takes form 

id \ I ~ ( Hdn{u)t,k) h(cn{ujt,k)-isn{u)t,k))\ ( <I'i(t) | ^^-j 

\'I'2(t) / \h(cn{ujt,k)+isn(u)t,k)) -Hdn{uit,k) J \^2{t) J 
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2 Independence of resonance on modulus k 

Let us turn to a dimensionless independent variable r = ut and substitute 
the dependent variables: 



^l(r) ] = [ / ] ( Mr) 



(6) 



where 



1 + CUT . . , , / 1 — cnr 



f = V CUT — isnr = y isign (snr) y . (7) 

The system in (4) takes on a different form through change in (6) 

,0jMr)\ ^(^^nr A \ ( Mr)\ ^ (g) 

V^2{r) J ^ \ h -^dnr J J 

in which detuning A is 

A = //-|. (9) 
The system in (8) may become similar to one studied by Shirley [51 through 

f V2 V 

a time independent orthogonal transformation ^ 

\ _ v/2 V2 
\ 2 2 / 

Let us note the mechanical-geometrical analog of the system in (8) while 
presenting the unknown functions as v^i(t) = x + iy, (/?2(t) = tt + zv, where 
X, y, u, V are real functions that are defined at the initial time and satisfy 
a non-linear system of differential equations with constant coefficients: 

x'^ + y'^ + + = 1, (10) 



h 

vdrX — udrV + ydrU — xdrV = — , (11) 

CO 



{drxf + {dry)' + {d^uf + (drvY + ^sn\T, k) = (12) 
ydrX — xdry + udrV — vdrU = — dn{T, k), (13) 



which describe the motion of a 4- vector {x, y, u, v) along the spherical surface 
(10) under two conservation restrictions (11), (12). 

The equation which is required to determine function (p2iT~) could be 
found from the system (8) 



drrM^) + { i—k'^snrnr -k'^sn^r + ) ^2{t) = 0, (14) 



where 

fi2j = /i2 + A2, (15) 

and is presented as a generalized Lame equation in the Jacobi form. 
Eq. (14) shows that the real u and the imaginary v components of function 
^2{'t) define a gyroscopic system with a parametric excitation of its intrinsic 
frequency ^ [Hj. 

Let /i(t) and /2(t) become a fundamental system that defines a general 
solution to Eq. (14) 

^2(r)=A/i(T)+5/2(T), (16) 

where A, B are arbitrary constants. The solution of the Cauchy problem 
(1) could be expressed through functions /i(t) and /2(t) using the initial 

conditions ( ) as 

V*2(0)/ 



where 



idrfi,2 + tfi,2dnT f,{0) F2 

F,,2 ^ ,ai - -^,a2 - 



T=0, 



^1 = -^^y^, b2 = ^|t=0, -d = -i{fldrf2 - f2drh)\r=0- 

a a uj 



:i8) 



Assuming that the wave function ( *^ ] is equal to ( ^ ) at the initial time, 

V*2(0)/ \oJ 

the solution for the wave function (17) at time r is given by 



] = ( f ° I I "l^l+^1^2 



(19) 



The probability of a transition requiring a spin flip over time r is 

Pi^_i{r,A,k) = \a,h + bj2\'. (20) 

Thus, the formula describing the transition probability could be expressed 
through functions /i(t), /2(t) using the formulae in (18) and (20) as 



Pi^_i(r, A, A;) 



Cj2 



/l(r)/2(0)-/2(T)/i(0) 



ifiir)drf2ir) - Mr)drMr))\,=o 



(21) 



The Rabi result P provided by Eq. (14) at A; = stipulates that /i(t) = 
cos /2(t) = sin and the transition probability is 

P._i(r,A,fc = 0) = ^sin^^r. (22) 

Eq.(14) is simplified at < A; < 1 in case of a sharp fundamental resonance 
with A = 0: 

drrM-r) + -^Mr) = 0. (23) 

Therefore, /i(r, A = 0, A;) = cos ^r, /2(r, A = 0, A;) = sin and Eq. (19) is 
solved explicitly 



^l{r) \ _ I feasor 



*2(r) / \ -if* sin 



(24) 



Obviously, the transition probability is independent of the k modulus and 
given by 

h 

Pi^_i (r, A = 0,0 < A; < 1) = sm^-T. (25) 

2 2 

The fundamental resonance is stable at any value of the k modulus with 
respect to consistent variations of the longitudinal and traversal magnetic 
fields. In other words, a distortion in the traversal field is fully compensated 
by a corresponding distortion in the longitudinal field . 

Knowing the wave function in (24) makes it possible to find the polariza- 
tion vector P defined by 

P, = (^(t)ai^(f)) {t = x,y,z). (26) 

A simple calculation produces 

P = {sn-frnHotsin-fmhot, -cnjrnHotsinjrnhot, cos-fmhot). (27) 
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The polarization vector satisfies the Bloch equation 

dtP = Jm[H,P], lm=^-Y- (28) 

3 Reduction to Heun equation 

A general case where both the detuning (9) and the modulus of the ellip- 
tical functions are different from zero requires changing the variable in Eq. 
(14) and taking advantage of the doubly periodicity of the elliptical func- 
tions. 

It is known that 

snir + iK') = , ^ cnir + iK') = -i f^^ , , < k < 1, (29) 
ksriT ksnr 

where the full elliptical integral of the first kind K' is 

/■f 1 

K' = = dip, k'^ + k'^ = 1. (30) 

° Y 1 — k''^si'n?ip 

Let us switch to a half- variable in functions snr , dnr using 



, 2sn^cn^dn^ , dn^^ — k'^sn^'^cv?^ , , 

snr = ^ ^ — r^, dnr = ^ h 31 

l-Fsn^^ l-Fsn^^ 

and introducing a new independent variable 

„t' , 1 (1 + k)sn^ - icn^dn^ 

z = sri^— = k-^- — ^ — — ^ ^, (32 

2 1 + ksn^^ ' ^ ^ 

we arrive at the algebraic form of an equation for function (P2{t) = y{z) after 
simple transformations: 

OzzV + { \ 7 H r 9zy+ 



z 



z-l 2; - Tz 



Bi B2 S3 , C*! , C'2 C3 \ ^ ^ 

+ T^-^. + 7- — TT? + — + ^-r + - — r h/ = 0' 



z 



2 



{z-lf {z-^f z z-l' z-^ 
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where 



Ai = A2 = As = ^, Bi = B2 = a-b, B3 = -(a + b), 

A , A2 

a = — , 



(34) 



n2 

Ci = 2(a - 6) - 2bk^ + (35) 



C2 = -2(a - 6) - (a + 6) + 7:^[«(^' - 1) " K^' + 1) + (36) 

/c 1 to 

C, = k\a + b)- T^Hk' - 1) - Kk' + !) + %■ (37) 

The values of Ci, C2, C3 are hnked through 

Ci + C2 + C3 = 0. (38) 

The functional form of the coefficients in Eq. (33) under the condition in 

(38) is necessary and sufficient to classify this equation as Fuchs' 

Eq. (33) may be simplified through substituting the dependent variable [31^1 

viz) = w{z)v{z), w{z) = z^{z - mz - ^y, (39) 

where p, q, r are both solutions to the indicial equations and characteristic 
exponents in the vicinity of the following points 

z = ai = 0, z = a2 = I, ^ = «3 = (40) 

The solutions may be found from the indicial equations: 

pip-l) + ^P + Bi = O^Pi = ^±\^-^\, (41) 

g(g - 1) + + = ^ g± = i ± 1^ - ^1, (42) 

s 1 1 , A 1, 

r(r-l) + -r + 53 = 0^r± = -±|— + -|. (43) 



The characteristic exponents in the vicinity of 2 = oo are determined from 
the indicial equation [7] 

3 3 

poo^poo _ 1) + (2 _ ^ + Y^iBk + = 0, (44) 

k=l k=l 

which takes the following form using (34-37), (40) 

- 1 ) + + 53 = ^ = 1 ± I A + 1 1 . (45) 

Thus, we obtain the Heun equation with real parameters [H] for function v{z) 
in the Klein-Bocher-Ince form P| 

where 

a = p++p + g + r, I3 = p'^+p + q + r. 
The accessory parameter qa is given by 

qa = ir + -p p-^- (48) 

It is easy to verify that the Fuchs' condition is satisfied: 

7 + 5 + e = a + /3 + l. (49) 

The Riemann symbol that characterizes the form of the solutions to Heun 
equation (46) is given by 

Q a z; qa \ ■ i^^) 

1-7 1-5 1-e (3 

It is worth noting that the substitution in (39) could be done eight ways 

pqr e {p+q+r+, p+q-r+, p+q+r^, p+q-V^, 

, (51) 

P-q+r+, p^q^r^, p^q^r^, p^q^r^] , 
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which sometimes allows to pick suitable characteristic exponents and ac- 
cessory parameter for specific calculations. Thus, the linearly independent 
solutions to Eq. (14) are written as fi = wvi, f2 = WV2, where vi, V2 belong 
to the fundamental system of solutions to the Heun equation (46). The prob- 
ability of the transition with a spin fiip is expressed through the solutions 
Vi, V2 as : 

w{T){v,{T)v2{0)-V2{r)vm) ' 

[w{T){v,{T)drV2{T)-V2{T)drV,{T))]\r=0 ' ^ ' 

In practice, solutions to the Heun equation are obtained through a re- 
duction of the parameter space for 7, 5, e, a, /3, qa taking into account that 
there are additional relationships between the parameters [9] besides the con- 
dition in (49). An analysis of parametric resonances induced by a non-zero 
detuning and distortion of a magnetic field will be presented elsewhere. 



A_.(r,A,^) = - 



4 Spin > \ 

Generalization for higher spin values is not difficult because the matrix 
in (17) is unitary due to the hermitian property of the Hamiltonian in Eq. 

(1): 

f f{a,F,+b,F,) f(a,F,+b,F,)\ =Jjl(^^^ (53) 
V/*(ai/l+bi/2) /* (12/1 +62/2)/ 

where Di{ip, 6, i/j) is a Wigner matrix. The Euler angles cp, 9, ip are 
determined from equations 

6 i 

f{aiFi + biF2) = cos-exp-{ip + ^), (54) 

6 -i 

/*(ai/i + bi/2) = i sin - exp— {Lp-ip), (55) 

Q 

sin^- = Pi^_i{T,A,k) (56) 

and are independent of the magnitude of the angular momentum J [TT] . 
Therefore, the transition probability over time r from the state with a pro- 
jection of angular momentum m into one with a projection m' for a particle 
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with a spin J is given by 

P„^„..(t,A,/c) = |-Dil,(jJ, 9, = |(J + m)!(J-m)!(J + m')!(J-"')!] 



4J^ 
COS - 

2 



(tan if^"^+^' 



v\[u — m + m')\{J + m — u)\(J — m' — v) 



(57) 



The Wigner matrix in formula (57) is 

D^mAv^^ V) =«"'""e'('"^+™'^)[(J + m)!(J-m)!(J + m')!(J-m')!]^ 

(sin 1)2'^-'^+'"' (cos |)2J-2^+'^-'»' 



— m + m')\{J + m — z/)!( J — m' — z/)! ' 

(58) 



and m, m' = {—J, —J+1, J— 1, J). 



5 Conclusion 

The study presented in this work may be useful in analyzing the results 
of interference experiments |il2j, improving magnetic spectrometers, and the 
field of quantum computing. 
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